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§ ; ABSTRACT 
(N I 

Q^, By radiation transfer models we show that the optical properties of grains 

! are poorly constrained by observations of reflection nebulae. The interstellar 

CN ! medium is known to be hierarchically clumped from a variety of observations 

^ ! (molecules, H I, far- infrared) . We have performed radiative transfer through 

^ ! four-tiered hierarchically clumped dust in a sphere surrounding a central star. 

! Our models have realistic power spectra of the projected density distributions 

^ ! (index ~ —3). The input parameters are the albedo (a) and phase parameter 

O ! (g) of the dust, the radial optical depth of the sphere averaged over all directions 

^ ! (tq), and the detailed random distribution of the dust clumps within the sphere. 

^ ! The outputs are the stellar extinction, optical depth, and flux of scattered light 

as seen from various viewing angles. Observations of a reflection nebula provide 
O I the extinction and scattered flux as viewed from one particular direction. 

I Hierarchical geometry has a large effect on the flux of scattered light emerging 

I from a nebula for a particular extinction of the exciting star. There is a very 

.!_H I large spread in both scattered fluxes and stellar extinctions for any distribution 

^ I of dust. Consequently, an observed (rext,Tsca) can be fitted by a wide range of 

_c3 I albedos. 

There are lower limits on a set by the scattered flux. As an example, in the 
best observed reflection nebula, NGC7023, a(1300A) must be larger than ~0.5 
if the scattered flux from Witt et al (1993) and a reasonable value for the optical 
depth within the nebula are adopted. However, the same observations can be 
fitted with a = 0.8 and 0.6 < g < 0.85, the entire range we considered. 

With hierarchical geometry it is not completely safe to determine even rela- 
tive optical constants from multiwavelength observations of the same reflection 
nebula. The problem is that the geometry effectively changes with wavelength 
as the opacity of the clumps varies. Limits on the implications of observing the 
same object in various wavelengths are discussed briefly. 

Henry (2002) uses a recipe to determine the scattered flux from a star with 
a given extinction. It is claimed to be independent of the geometry. It provides 
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considerably more scattering for given dust optical properties than our models, 
probably leading to an underestimate of the grain albedos from the UV Diffuse 
Galactic Light. 

1. Introduction 

Interstellar dust is poorly understood in spite of its importance within galaxies. The 
mean extinction law (i.e. scattering plus absorption) for the diffuse Galactic interstellar 
medium (ISM) is reasonably well known, but the individual absorption and scattering prop- 
erties are still controversial. Many other basic properties (chemical composition, particle size 
distribution, particle shape distribution, alignment of the grains in the magnetic field, and 
others) are also poorly known. 

The properties of dust for each scattering event arc characterized by two parameters: 
a(A), the albedo or fraction of the extinction that is scattering, and g{\), the mean cosine 
of the scattering angle, averaged over the phase function for single scattering. Geometrical 
uncertainties prevent the determination of higher moments of the scattering pattern beyond 
these two. The a and g are estimated from either reflection nebulae excited by a central 
star, radiation from a dusty globule reflecting Galactic light, or the "diffuse galactic light" 
(DGL), the general stellar radiation reflected by the dust distributed throughout the ISM. 

Almost always, nebulae have been interpreted with smooth (non-clumpy) models, some- 
times with a density gradient. Radiation transfer differs in smooth and clumpy media be- 
cause photons can escape between the clumps. Uniformity was a reasonable assumption, for 
three reasons, (a) It is the simplest, and Occam's Razor has great appeal, (b) Hierarchical 
clumping had not been convincingly demonstrated by several means, (c) The computational 
resources required for radiative transport through hierarchical clumps have only recently 
become widely available. 

Any estimate of the albedo is affected by the hierarchical clumping of the ISM. These 
clumps occur on scales from AU to tens or hundreds of parsecs, as observed in various 
molecules (CO, formaldehyde, and OH), ionic species, H I, and 100/im (Crovisier & Dickey 
1983; Green 1993; Moore & Marscher 1995; Vogelaar & Wakker 1994; Elmegreen & Falgarone 
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1996; Abergel et al. 1996; Elmegreen 1997; Heiles 1997; Heithausen et al. 1998; Falgarone et 
al. 1998; Chappell & Scalo 2001; Welty & Fitzpatrick 2001; Faison & Goss 2001; Andrews, 
Meyer, & Lauroesch 2001). Other galaxies (Stanimirovic et al. 1999 for the SMC; Westpfahl 
et al. 1999 for the M81 group) show the same phenomena in H I. The best observed reflection 
nebulae are clumpy (e.g. Sellgren, Werner, & Dinerstein 1992; Martini, Sellgren, & DePoy 
1999; Knauth et al. 2001). We have examined the images in the SIMBAD database of the 
reflection nebulae with which we are acquainted (about a dozen); all appear to be clumpy. 
The observed thermal pressure in the ISM as judged from C I (Jenkins & Tripp 2001) varies 
by over an order of magnitude within a sample of stars against which the column densities of 
C I hnes can be determined. The hierarchical structure of the ISM is a natural consequence 
of tiirbnlcncc, which is scale-free. Turbulent models of the ISM generate hierarchical density 
structure (e.g., Norman & Ferrera 1996). 

Boisse (1990) determined the transport of radiation through a two-phased clumpy 
medium with isotropic scattering. Witt & Gordon (1996) extended the calculations to 
anisotropic scattering in a centrally illuminated sphere. They investigated the effects of 
changing various parameters such as the contrast between the density contrast between the 
phases or the sizes of the clumps relative to the radius of the sphere. Many of the qualitative 
effects we find from varying parameters for hierarchical models are discussed by them as 
well. As the volume fraction of the dense phase was increased, their models varied from 
widely separated dense clumps up to continuous structures with small holes. Since all cells 
have either of two densities, we will refer to this type of clumping as "two-phase models" . 
By contrast, our hierarchical models have an almost continuous distribution of densities. 

Witt & Gordon (2000) have investigated radiative transfer in galaxies by means of a 
specific two-phased ISM, one being 100 times as dense as the other and occupying 15% of the 
volume. In order to compare two-phase models with hierarchical, we adopt this recipe, with 
16 cells along the radius of the sphere (Witt & Gordon used 15). All of our two-phase models 
assumed (a, g) = (0.6,0.6), which arc values typical for optical wavelengths suggested by 
observations (see plots in Witt & Gordon 2000). 

The fact that the ISM is hierarchically clumped is important for the propagation of 
scattered light, since such a structure has relatively open spaces through which the radiation 
can move rather freely. The main thrust of this paper is not to try to put limits on the grain 
properties of real reflection nebulae, but to show the variation among hierarchical models 
as they are viewed from various angles. We will not be very concerned with the averaged 
properties of the models, since each real object is viewed from only one direction. 
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2. Hierarchically Clumped Models 

The density structure of the ISM is often described as "fractal", meaning self-similar 
on all size scales. We will consider models that are hierarchically clumped instead, meaning 
that they are self-similar over a limited range (about a factor of ten) in sizes. We use a 
procedure similar to that in Elmegreen (1997): (a) Consider a "supercube", a portion of 
which will represent a spherical reflection nebula. The supercube of size L on a side consists 
of 64 cubical cells stacked along each dimension. For each cell we determine the local density 
of dust, as explained below, (b) Place N points randomly within the supercube. We used 
= 32. (c) Randomly cast another points, all within a distance L/(2A) in each Cartesian 
axis from each of the points cast in the preceding round. Here the distance A is related to 
the "fractal dimension", D,hj D = \og{N)/ log(A). Allow any of the points that fall outside 
of the supercube to remain there, (d) Repeat procedure (c), above, twice more, so that the 
total number of points cast is N'^. (e) Shift the points outside of the supercube to within it by 
translating each Cartesian coordinate outside of the supercube by L until it lies within. This 
procedure reflects the points so that if they originally fall off the left side of the supercube 
they reappear the same distance from the right within the supercube, and correspondingly 
for up/down and front/back. The density within the supercube is then proportional to the 
number of points within each cell, (f) Inscribe a sphere within the supercube and place a 
point source of radiation at its center. The constant of proportionality between the number 
of points in each cell and the optical depth within the cell is chosen to make the radial optical 
depth, averaged over all directions, be tq. 

The main parameter of this procedure is D, which describes the degree of clumping. We 
considered D = 2.3 and 2.6, in the range observed in actual clouds. Another observable pa- 
rameter is /?, the exponent of the power of the projection of the density within the supercube 
onto the plane of the sky^, so that P{k) oc k~f^, where k is the wavenumber. Observations 
(e.g. Stanimirovic et al. 1999) show /3 ~ 3. With our recipe we find /3 = 2.8 with D — 2.6 
and P ~ 2.3 with D — 2.3. However, none of the conclusions in this paper depend upon 
which set of hierarchical models we consider. Many different models can have the same p. 
No one parameter, either D or /?, can completely describe even the projection of the density 
distribution of the ISM onto the sky. 

The above procedure with the parameters we adopted results in clouds that are quite 
full of holes, with ~15% of the projected density distribution being almost blank. The 



"^Thc structure function F{5r) is defined by F{Sr) = (area)""^ / / /(r)/(r + 5r)dxdy, and P{k) is its 
Fourier transform, with k = |<5r|~^. Here r and 5r are vectors, and {x, y) Cartesian coordinates, in the plane 
of the sky. 
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large-scale ISM as shown by H I is similar if a low column density of H I (< 5 x 10^^ H 

atoms cm~^) is taken to be "blank" (i.e., a column density much smaller than the average). 
Elmegreen (1997) has discussed the emptiness of the ISM, with the large-scale structure 
in mind. Images of actual reflection nebulae (see Sellgren, Werner, & Dinerstein 1992 for 
NGC 2023 and NGC 7023) show strong filaments and structure, but there is material at all 
points within the nebula. For this reason, we sometimes add a uniform density to all cells 
in the supercube. This constant density does not change the value of since the Fourier 
transform of a constant is a Dirac 5-function at the origin, k — 0. Choosing a constant 
density provides a minimum projected density. 

The detailed placement of the clumps of dust also infiuences the results of our radiative 
transfer calculations. The clumping is mainly formed by the first round of random casting 
of 32 points that are subsequently spread by three more rounds of casting points in their 
vicinities. The locations of the first round of points is uniquely determined by the integer 
seed of the random number generator^, and the locations of all subsequent points follow from 
the initial value of this integer in a complicated but unique way. The response of the model 
nebula to the central star is strongly affected by the precise placement of the dust relative 
to the star and, thereby, to the value of the initial seed. We considered suites of models that 
differed only in their initial seeds. Our procedure amounts to assuming that the star does 
not affect the density distribution of the dust in its immediate vicinity, so the initial seed is 
allowed to control the placement of the dust both near and far from the star. To test the 
sensitivity of results to this assumption, we also considered models with a cavity in the dust 
distribution within the inner 10% of the radius of the sphere. 

The radiative transfer was performed by the Monte Carlo code described by Wood & 
Reynolds (1999). It involves considering photon packages propagated in a random direction 
from the central star. We assumed a and g, along with the scattering phase function of 
Henyey & Greenstein (1941). In most cases we used 5 x 10^ photons for each model, after 
checking for a few cases that the results were the same as from a model using 2x10^ photons. 
Unlike the situation for placing the dust, the initial seed for the radiative transfer makes no 
difference, as expected from such a large number of stellar photons propagated in random 
directions. 

We calculated the fiuxes of scattered radiation and also of starlight, both of which 
depend on the angle from which the sphere is viewed because of the strong randomly placed 
clumping. Both fluxes are expressed relative to the stellar flux that would be observed from 
the unobscured star. We express the fluxes as extinction and scattering optical depths. Text 



^We used ran2(iseed) as described in Press et al. (1992). 
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= — ln(F*/F*o) and Tgca = — ln(Fsca/-^*o), where F*o is the flux the star would have if there 
were no reflection nebula. Each flux or optical depth depends on the direction from which the 
sphere is viewed. Wc considered 18 evenly spaced values of cos(^^) and 36 values of 0, where 9 
is the polar angle of the sphere and (j) the azimuthal. Fgca and can be determined directly 
(with difficulty!) in real nebulae. F^o can be determined if there is enough wavelength 
coverage to determine the extinction from extrapolation to large wavelengths, from colors if 
the reddening law is assumed, or if the total FIR flux from the nebula is measured. Otherwise, 
the ratio of the stellar to the nebular flux gives the difference of Text and Tgca- We now ask, 
how well can we determine (a, g) from an observation of (Text, ^sca)? 



3. Results 

As regards radiative transfer, hierarchical clumping introduces major differences from a 
uniform distribution or, to a lesser extent, from two-phase models. The main properties of 
hierarchical models are: 

(a) When viewed from various directions, hierarchically clumped nebulae exhibit Text 
and Tsca spanning a wide range. A real reflection nebula represents viewing a collection of 
dust from one particular direction. We can calculate mean properties of models as averages 
over all viewing directions, but mean properties are relevant only if a large collection of 
nebulae (such as in another galaxy) are observed. Angle-averaged properties of any model 
have very limited applicability for interpreting a given reflection nebula. The two-phase 
models models of Witt & Gordon (1996) also show a wide range of Text at a given tq. 

(b) As expected, the global parameters of the models, such as Tq and D, have a significant 
influence on the radiation transfer. We considered both purely hierarchical models and some 
with 33% of the mass in a constant density sphere upon which the clumps are superimposed. 
This uniform background density can have an important effect on the radiation transfer in 
optically thick nebulae, when the diffusion depends more upon the transparent regions than 
on the opaque. This rule is familiar in connection with the use of the Rosseland mean opacity 
in the optically thick limit of radiation diffusion in stars. Whether such a uniform layer exists 
in reflection nebulae is not clear. The regions between clumps on large (molecular cloud) 
scales have a vastly lower density than the mean (e.g., Elmegreen 1997). 

(c) Models with the same global parameters (tq, -D, a, g), but differing in the placement 
of the dense clumps relative to the central star, show strong variations in even the aver- 
aged fluxes. There are signiflcant differences between hierarchical models in which the star 
happens to lie within a dense clump of dust, in contrast to those in which it falls within a 
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void (or if there is a central hole). If the star happens to lie within a void, there are lightly 
reddened paths that reach from the star to the edge of the sphere. If the star is almost 
unattenuated as seen from a particular direction, the nebula would be interpreted as having 
the star in front of the dust, except that there can be a significant amount of scattered light 
that would not be there if the star were truly foreground. In these cases, the albedo derived 
from uniform models would be > 1. 

The power law index, f3, of the projected density distribution of our two-phase models 
is ~0, reflecting the fact that their projected density is non-hierarchical. The two-phase 
models always provide a smaller range in Tscaid-ct)) (about ±0.12) than the hierarchical, but 
their spread in Text is comparable to hierarchical models at low tq (< 2). The spread at large 
To decreases because the interclump medium becomes opaque, so that in the ultraviolet (UV) 
two-phase models make a prediction of a that is rather independent of viewing angle. 

Figure 1 shows Tgca plotted against Text, for the case tq = 2, D = 2.6, a = 0.6, g = 0.6 
for a hierarchical distribution of clumps with an initial seed for the density distribution that 
provides a typical (defined below) placement with respect to the star. Each point represents 
a particular viewing direction. One sees a large spread in the values of Text- The values for 
Text for this model range up to 6.8, and Tgca to 2.4. Thus, this one model of clumpy dust can 
produce a wide range of stellar fiuxes and scattered light. 

We represent the direction-averaged escaping fiux (not mean optical depth) by (Text), 
defined by 



where iV^ ^ is the number of bins we are using (= 18 X 36), all of equal solid angle. If the star 
is embedded within a clump, the uniform dust component is of reduced importance because 
scattering occurs within the clump. Similarly, (rgca) represents the angle-averaged scattered 
radiation, with Tgca replacing Text in the above relation. The filled square in Figure 1 shows 
(Text) and (Tsca), which are useful for contrasting the results of models with the same value of 
To but different initial seeds. As the figure shows, the ((Text), (Tgca)) point is of very limited 
value in guessing what actual optical depths a given nebula might show if it were viewed 
from various angles. In this model (Text), which describes the average optical depth of the 
stellar flux that escapes from the nebula, is only O.Stq, the averaged optical depth through 
both the clumps and interclump material. It is sobering to recall that only a handful of 
reflection nebulae have been well observed, while each point in the figure represents a view 
of a single model with a particular clumping and a, g, and tq. 

The solid lines in Figure 1 are the tracks of a spatially uniform dust distribution for 
the same g {— 0.6) and various values of a, as labeled. These lines were computed with 




(1) 
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the surprisingly accurate analytical approximation (the scattered flux is underestimated by 
<3% for Text < 4) given by Code (1973). Most reflection nebulae have been interpreted using 
these uniform models, in which case Tgxt is the true optical depth of the nebula. For our 
hierarchical model, in about 15% of viewing directions we would conclude that the star is 
in front of most of the dust because the uniform models would give a > 1. In 40% of the 
directions we would conclude that a > 0.6, the value assumed in the hierarchical model. 

Two UV observations of NGC 7023 arc shown in the figure. The central star, HD 200775, 
has foreground extinction from ordinary ISM in addition to the extinction provided by the 
nebula (see Witt ct al. 1993 for discussion and other references to observations of NGC 7023). 
The circled asterisk marks the nebular parameters derived from Voyager-2 and Hopkins 
Ultraviolet Telescope (HUT) data if 40% of the total extinction arises within the nebula, 
and the circled dot if the nebula provides 50%. The dashed fines will be discussed in §4 in 
connection with the UV DGL. 

Figure 2 is the same model as Figure 1, except 33% of the dust is in uniform distribution 
upon which the hierarchical clumping is superimposed. The axes of the figures are the same 
as Figure 1 to emphasize that a uniform density increases the scattered light, so long as its 
optical depth along a radius is not too large (< 1), and imposes a minimum value of Text- 
The range of Tsca for a given Text is smaller than in Figure 1, but this is an artifact of the 
uniformity of the extra component that is the only difference in the two figures. 

Wide areas of the (Text, Tgca) plane are covered by each hierarchical model, so a given 
(^ext, "^sca) observation can be covered by a wide variety of models. By contrast, uniform 
models have very limited dependence on g, and so contain only one effective parameter, a. 
The ambiguity of hierarchical models is suggested by Figure 3. It has the same axes as 
Figures 1 and 2, but only outlines around the areas covered by the points of five hierarchical 
models. The solid line has a high a (= 0.8), moderate g (0.6), and low Tq (0.5). The maximum 
values of Text ai'c comparatively small because of the low tq. The dashed line encloses the 
points for the same model with 33% of the dust mass in a constant density; the scattered 
points are lower (there is more scattering) because there are no almost empty paths through 
the nebula. The dot-dashed and dotted fines are the boundaries of the points in Figures 1 
and 2, which arc for (a, g^ D, tq) = (0.6, 0.6, 2.6, 2). The long dashed lines enclose the very 
wide boundaries for (a, g, D, tq) = (0.8, 0.85, 2.3, 4). We can have large g or moderate, 
large a or moderate, and large tq or sm,all and still fit the observations. However, we will see 
that a minimum a of ~0.5 is required to produce enough scattered light. 

Figure 4 shows the averaged ((Text), {Tsca}) values for all of the 21 different initial seeds 
that wc tried. Figure 4a shows purely hierarchical models; 46, models with 33% of the dust 
in a constant distribution. The parameters are the same as in Figures 1 and 2: a = g = 
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0.6, To = 2. The NGC7023 points are as before. The open squares are models with D — 
2.3 instead of 2.6. Within each panel, the D — 2.3 models are higher because the radiation, 
both scattered and stellar, can escape more easily from the more strongly clumped structure. 
The lines show uniform models, with albedos marked in 4a. The uniform model appropriate 
to the actual albedo assumed in the models (a = 0.6) is just below the bottom of the figure. 
The dashed line is the uniform model with a = 0.5 and g = 0.4 instead of 0.6. The differences 
are not large in comparison to the effects of the other parameters. The effects on hierarchical 
models of changing g are similar. 

Perhaps the most striking difference between the two panels is the lower values of (rgca) 
for hierarchical models with 33% of the dust in a uniform component. The increase of 
scattering from the dust between the clumps causes the decrease in (rgca) and greatly reduces 
the differences between models with different spatial distributions of dust. In either panel, 
the points from various viewing angles of the two values of D are completely intertwined. 
Reflection nebulae are poor diagnostics of D as well as other properties of the ISM. 

The three open circles in Figure 4 are two-phase models with three values of the initial 
seed. Two of these values show the extrema in ((Text), (Tgca)) among the 21 initial seeds that 
we tested. The tightness of the mean optical depths of the two-phase models models shows 
the importance of hierarchical geometry, as opposed to simple two-phase clumps. The (Tgca) 
of the two-phase models models are similar to those of the hierarchical models with uniform 
dust. The contrast of both with purely hierarchical models illustrates the importance of 
dust-free regions (in real space, possibly caused by extensions of hot, low-density material 
into the nebulae). 

Figure 4 shows what we meant by saying that Figure 1 was produced by a typical 
hierarchical model - one with a typical (rcxt)- Those with large (rcxt) have the star embedded 
within dusty material, and a relatively low central dust density leads to a low (Text)- 

Differences among models, and errors arising from interpreting reflection nebulae with 
uniform models, increase with optical depth. Figure 5 shows the same as Figure 4 (a = 
0.6; g = 0.6), except that Tq= 4 instead of 2. The open squares are for various initial dust 
distributions with D =2.3, with 33% of the density uniform; the filled squares, the same 
with D — 2.6; open triangles, purely hierarchical density, D—2.3; filled triangles, purely 
hierarchical with D — 2.6. The circles enclose the entire range of our two-phase models 
models. 

Figure 5 shows an important result: optically thick hierarchical models can show very 
modest values of (Text), but interpreting the nebula with uniform models will likely result in 
a severe underestimate of the albedo unless there is material between the clumps to scatter 
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light efficiently. For the purely hierarchical models with D — 2.3 (open triangles), there is 
very little scattered light between the clumps (i.e., large Tgca) and limited stellar extinction, 
resulting in albedos predicted from uniform models that are as low as 0.17. The spread 
among the two-phase models (circles) is relatively small because at large optical depth they 
become statistically uniform. Since they are clumpy, the uniform-model albedos they suggest 
are significantly lower than the 0.6 used in the actual two-phase models. 

The squares in Figure 5 have 33% of the dust mass in a constant density, so that the 
minimum optical depth is 4/3. This addition greatly increases the scattered flux and makes 
the difference between D = 2.3 and 2.6 almost unimportant. The albedo that would be 
derived from uniform models is increased to ~0.5 over the purely hierarchical models, still 
short of the correct value of 0.6. 

The spread of the (rcxt) and (rgca) in Figures 4 and 5 show that the spcciflcation of 
To and D (along with grain properties) does not specify the scattered and stellar fluxes 
from reflection nebulae in even a statistical sense. By contrast, D constrains the statistical 
properties of the projected brightness distribution for both line and continuum emission 
from the gas. The large variations in properties of reflection nebulae with viewing angle 
come about because the radiation arises from the geometry of the statistical distribution 
of the ISM relative to a point source rather than to itself. Emissions form gas involve an 
integration over larger regions of more equal weight. 

It is unrealistic to suppose that reflection nebulae are spherical with a central exciting 
star. We have considered only this geometry because we feel that it provides rather general 
results for strongly forward-throwing scattering (say, g > 0.5). In this case, the dust on 
the near side of the nebula is likely to provide most of the scattered light. We consider our 
models as representing only this near-side dust. 

4. Discussion 

4.1. Reflection Nebulae 

Our hierarchical models, with their overall spherical distribution of clumps, are primarily 
aimed at interpreting reflected light from reflection nebulae illuminated by an exciting star. 
We will discuss reflection of the diffuse Galactic radiation held from dusty globules at the 
end of this section. 

We now examine the limits on a imposed by the amount of scattered light observed in 
reflection nebulae. Since real nebulae show the illuminated cloud as it is viewed from only one 
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direction, we must consider the spread of (Text, ^sca) from any particular dust distribution. 
We seek Omax and amin, the minimum and maximum values of a that can fit the observations. 

For illustrative purposes, we consider the 1300 A observations of NGC 7023 with 50% of 
the stellar extinction assumed to be within the nebula (Witt et al. 1993; the circled dot in 
Figures 1 and 2). In Figure 1 it lies in the middle of the envelope of the (Text, Tgca) points from 
various viewing angles. If we decreased the a of the model, the scattered flux at each direction 
would become fainter and Tgca larger. The stellar fluxes would be unchanged because they 
depend only on extinction, not scattering. Thus, the pattern of points would move up in 
the diagram, and the point could no longer be fitted. To find the amin for NGC 7023, we 
must vary tq and initial seeds (dust distributions) and find the minimum a for which the 
envelope of individual points lies just above the observation. A similar procedure leads to 
Qmax- Of course, it is improbable that the actual o is as low as Omin 

because the observed 

(^ext, Tsca) would probably not lie at the edge of the (Text, Tgca) distribution as seen from a 
random direction. 

We have performed this exercise and found the following: 

(a) All of our models required an albedo of > 0.5 to fit the lower point of NGC 7023. 
Varying parameters and geometries to obtain the absolute minimum a (0.46; see below) 
produced only changes in details. This result is robust because the observed Text is near the 
most efficient value for producing scattered light, as suggested by the uniform models (see 
the minima of Tgca in the curves in Figure 1). 

(b) The scattered light is near its maximum when tq ~2, where the mean Text is about 
unity, so Omin is best estimated for our hierarchical models near this value. The maximum 
scattered hght for uniform models occurs at Tq ~ 1.1 - 1.4, increasing slowly with the albedo. 

(c) The smallest amin we found for the circled dot in Figure 1 is 0.48. It occurs for tq 
~2, D — 2.6, and a particular seed whose ((Text), (Tsca)) point is low in Figure 4. The effects 
of g are minimal on the locus of the lower envelope of the individual points; the minimum a 
is achieved for g > 0.5. We tried such a wide range of models that we doubt that there is 
an amin drastically lower. 

(d) No useful estimate of amax is possible because inefficient scattering can be produced 
by extreme clumping (as small a D as allowed) and with no uniform component. For D — 
2.3, an a ~ 1 is needed to fit NGC 7023 for some distributions of clumps. For D — 2.6, amax 
is about 0.8. The dust geometry that produced amin did not allow a to be above 0.6, but it 
was chosen to have especially efficient scattering. Alternatively, very optically thick models 
that have low stellar extinctions as seen from some viewing angles can fit observed fiuxes 
with very high albedos. For instance, we have fitted the NGC 7023 point in Figure 1 with 
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two-phase models, tq — 8, a — 0.85. We have no doubt that many optically thick hierarchical 
models could fit the observations as well. Perhaps such extreme models could be ruled out 
with far-infrared (FIR) fluxes, although the low absorption tends to compensate for the high 
optical depth. Furthermore, the FIR arises from dust surrounding the star in all directions, 
while the scattering is mainly from dust close to the line of sight. In any case, these models 
are hardly reasonable on physical grounds. 

Voyager 2 and HUT observed NGC 7023 at lOOOA (Witt et al. 1993). Uniform models 
predict a ~0.45 if 40% of the stellar extinction occurs within the nebula. Hierarchical models 
can fit the observation (not plotted, but in the upper right corner of Figures 1-3) with a 
— 0.8 with To > 4. Of course, we make no claim that such a high a is correct. On the other 
hand, we can fit the point with a — 0.4, lower than the uniform models. Once again we see 
that the observations are subject to a very ambiguous interpretation. 

Burgh, McCandliss, & Feldman (2002) used a uniform model for NGC 2023 with Text 
= 2.4. They derived an albedo of O.SQtom at 1345 A. Our hierarchical model shown in 
Figure 1, with a — 0.6, could accommodate their observations, and others that scatter hght 
inefficiently could have a ~ 1. 

Figure 6 shows the albedos derived from uniform models, using the scattering and 
extinctions from the hierarchical model (a, g, D, tq) = (0.6,0.6,2.3, 1), plotted against the 
stellar extinctions at each viewing angle. The albedo used to generate the scattered fluxes and 
extinctions, 0.6, is marked by the dashed line. The distribution of points for the hierarchical 
model with tq = 2 is virtually identical. We see that the minimum a derived from uniform 
models is '--'0.32; the average for 0.6 < Text < 1-6 is '--'0.42. At moderate extinction, leading 
to bright scattered flux, uniform models underestimate a. At small Text uniform models 
overestimate a. A hierarchical model in which the star happens to be within a clump has 
smaller errors. Adding a constant density decreases the errors at Text ~ 1 so that the average 
a is ~ 0.5. Of course, there is no error in the unrealistic case of the density being uniform 
throughout. 

An important question is, how reliably can we determine a(A) from observations of a 
particular nebula covering a significant range in wavelengths, implying differing opacities or 
optical depths? In this case, we might hope to find relative optical properties of the grains 
accurately even with uniform models. Figure 7 shows the results for two wavelengths for 
which the extinctions (both tq and Text) differ by a factor of 2: for example, (1000 A) / (1540 A) 
if the extinction law is the Galactic average (Cardelli, Clayton, & Mathis 1989) and Ry — 
A{V)/[A{B) - A{V)] = 3.1. The figure displays the ratio of a(To=2) to a(To=l), both 
albedos derived from uniform models, for the same hierarchical model as shown in Figure 
1. The abscissa in the figure is Text for the lower extinction. All points in the figure were 
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computed with the a — 0.6, so the "true" value of the albedo ratio is unity. We see that 
uniform models have difficulty in predicting even relative values of a(A), since the spread 
in a(ro = 2)/a(ro = 1) is ~40%. The figure gives the impression that if r^xt is > 1.5, 
uniform models would derive a greater albedo at short wavelengths (say, 1000 A) than at 
longer (1540 A for a factor of 2 in the extinction). For NGC 7023 the opposite is found to 
be the case (Murthy et al. 1993; Witt et al. 1993). From Figure 7 we would conclude 
that a(1000A)< a(1540A) is robust, since Witt et al. (1993) suggest Text(1540A) ~(1.5 - 
2). However, models with the star embedded in a clump have the albedo ratios shift from 
mostly < 1 to mostly > 1 at Text(1540A) ~ 2.5, considerably greater than shown in the 
figure. These models would produce an illusory result a(lOOOA) < a(1740A). 

The basic problem is that the importance of a clump depends upon its optical depth. 
If a clump becomes optically thin, its existence no longer very significant to the radiative 
transport. Similarly, increasing the optical depth of an optically thick clump has a limited 
effect on the radiation. While the density structure within the ISM is independent of the 
wavelength, its geometry as regards radiation transfer is not. Since the clumping produces 
the large spread in Tsca.{0, 0), it is not surprising that there is a wide spread in the importance 
of the clumps among various viewing angles. 

The problems of determining the variation of albedo at various wavelengths are small 
if the ratio of optical depths is close to unity. Statements that the albedo is smaller at 
the 2175 A feature than on either side of it (e.g., Witt et al. 1992) seem robust, besides 
being completely plausible physically. At nearby wavelengths A{X)/N{}i), the extinction per 
H nucleus, can be the same at three wavelengths, and comparison of the relative albedos 
between those wavelengths is completely safe. 

We have seen that several properties of reflection nebulae change if the exciting star is 
embedded within a rather dense clump of dust. There is an increase in the scattered light; 
the lowest values of (rgca) for purely hierarchical models in Figure 4a (i.e., the large amounts 
of scattered light for the purely hierarchical models) belong to such cases. The central clump 
ensures that the angle- averaged extinction, (Text), is larger than for most models with the 
same tq (i.e., the points in Figure 4a tend to have large (Text) if (^sca) is small). Figure 46 
shows that the addition of a constant density component makes the effects of a central clump 
rather unimportant because the radiation scattered by the uniform dust dominates. Another 
effect of a star being embedded in a dense clump is that the errors of the albedo derived from 
observations of the same object at different wavelengths, illustrated in Figure 6, are reduced 
because more of the scattering arises from the central clump. For such a clump, the opacity, 
but not the geometry, changes with wavelength. Finally, there would tend to be more FIR 
radiation from an embedded star because the grains are relatively close to the star, so they 
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are absorbing a relatively large stellar flux, and the mean temperature of the grains would 
be relatively large for the same reason. 

Galactic radiation scattered by isolated globules seen away from the Galactic plane 
provides another diagnostic of grains (Mattila 1980). Witt, Oliveri, & Schild (1990) have 
analyzed such a globule and give references to other examples and discussion of the process. 

The density is concentrated towards the center of such globules. The reflected intensity rises 
outward from the center, peaks, and falls off to the sky brightness at the edge, showing that 
the scattering is forward-throwing. The reflected intensity can be related to the incident 
radiation predicted for both stars and the DGL out of the Galactic plane. Results were that 
at 4700 A grains are forward-throwing [g ~0.8) and very efficient at scattering (a ~ 0.8). 
The globule was modeled with a very centrally condensed density distribution that is smooth 
and spherically symmetric, as seemed completely appropriate at the time. Since turbulence 
driving clumps in the diffuse ISM may decay when driving into a steepening density gradient, 
perhaps the ISM within globules is relatively free from clumping. We cannot assess the 
uncertainties in the grain properties derived from globules until we know more about the 
clumping within them. The major uncertainty would then be the prediction of the radiation 
incident upon them from the hierarchically clumped Galaxy. 



4.2. The Diffuse Galactic Light 

There is no doubt about the applicability of hierarchical clumping to the DGL. Obser- 
vations such as those at 1740 A (Schiminovich et al. 2001) are based on large portions of 
the sky. We have seen that uniform spherical models of reflection nebulae underestimate a if 
applied to optically thick situations (large Tq, but not necessarily large Text)- Our reflection 
nebula models are not ideally suited for interpreting the DGL in the UV, since distributions 
of the dust and exciting stars are skewed towards the Galactic plane. However, our results 
have strong implications as regards the interpretation of the DGL as discussed in Henry 
(2002, hereafter H02). In that paper it was concluded that a < 0.1 for A < 1400 A. 

H02 performed radiative transport calculations by using B, V magnitudes and the spec- 
tral type for each star in the entire Hipparchos Input Catalog. The reddening, E{B — V), 
follows from {B — V) and spectral type. The extinction Ax is E{B — V) Ry (Ax/Ay), where 
Ry is usually taken to be 3.1. The Text for each star is A{X)/1.0S6. 

The heart of H02 is its unique treatment of radiative transfer. The scattering is taken 
to be 

exp(-Tsca) = a(l - e-^-*) exp{-(l - a)[Text + ln(l + e""-*) - ln2]} , (2) 
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independently of the geometry. However, Witt, Priedmann, & Sasseen (1997) have shown 
from clumpy models that the distribution of cloud optical depths can strongly affect the 
interstellar radiation field. 

At low optical depths, the H02 relation between Text and Tgca approaches the uniform 
model. The dotted hues in Figure 1 show the H02 models for three albedos identified by 
the nearby solid lines. We see that the H02 assumption predicts more scattered light than 
uniform models, and a lower albedo follows. 

For A< 2000 A, the Text for most hot stars is probably rather large, since the mean A{V) 
is ~1 magnitude kpc~^ and the ultraviolet is 3 - 5 times larger. For Text = 1-5, Figure 1 
shows that the H02 model with a = 0.6 predicts a scattered flux ~exp(0.4) = 1.5 times larger 
than that "typical" hierarchical model, chosen for illustrative purposes only. The conclusion 
of H02 that a(lOOOA) is very low ("perhaps 0.1") is perhaps understandable. Such a low 
albedo is not compatible with the observations of NGC 7023 by Witt et al. (1993), in which 
the scattered fiux was larger than the stellar, or Murthy et al. (1993), which suggests that 
a(1150 A) ~ 0.55 if 40% of the extinction is within the nebula, as assumed by Witt et al. 
(1992). 

Murthy, Henry, & Holberg (1991) used Voyager 2 to search for DGL for A < 1300 A. 
They failed to detect any, possibly confirming the H02 picture. Our refiection nebula results 
commonly show that the same nebula can show very low and rather high scattered fiuxes as 
seen from different angles. We strongly suspect that the same property apphes to the DGL. 
Our general result from hierarchical models is that photons, especially the unscattered ones 
coming directly from the star, can leak out of dust distributions far more easily than either 
uniform- density or H02 models would predict, so a dark sky does not necessarily imply a 
low albedo. 



5. Final Remarks and Summciry 

We have arrived at a pessimistic, but we feel realistic, assessment of the ability of 
refiection nebulae illuminated by embedded stars to serve as diagnostics of interstellar grain 
properties. It is better to realize what we do not know, rather than believe that we know 
more than we do! Even if the scattered light from a reflection nebula and the extinction 
suffered by its exciting star are observed perfectly, the geometry of the dust is vital in the 
interpretation of the observations. Hierarchical clumping, "clumps within clumps within 
clumps" , has a very large impact on the interpretation of both reflection nebulae and the 
DGL as compared to smooth density structures. Hierarchical structures, demonstrated by 
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CO, 100 //m, various other molecules, and H I observations, make any determination of the 
optical properties of grains from reflection nebulae very difficult. Our models had plausible 
power spectra for the autocorrelation function of the projected densities, with an exponent 
of the power of ~ —3, unlike either uniform models or those that arc clumped on a single 
scale (exponent ~0). The exact geometry of the dust is very difficult to determine because 
it depends on the placement of hierarchical clumps of dust around the star. 

We have made no interpretations of any reflection nebulae. It is now relatively easy 
to use Monte Carlo transfer codes to predict the scattering (and polarization) for objects 
with a known geometry. Inverting observed intensities to determine the correct geometry, 
though, is very difficult, even with observations with good spatial resolution taken in multiple 
wavelengths. Our hierarchical models have been schematic and do not begin to exhaust the 
possibilities of the real geometry of any reflection nebula. Observations in the infrared and 
line emissions with good spatial resolution may eventually put useful limits on the actual 
geometry of the gas and dust, but the problem of inverting the data will always be difficult. 

An interesting possibility is that isolated Bok globules may not be strongly clumped, 
so reflection of Galactic light from them can diagnose grain properties. The penetration 
of turbulence into them might be weaker than its percolation through larger structures. 
Nonclumpy globules can serve to diagnose grains because they scatter incident Galactic 
radiation (see Witt et al. 1990 for references and an explanation). 

Each of our hierarchical models shows a very wide range of extinction optical depths 
(Text) and scattering (rsca) when the nebula is viewed from various directions. The Text and 
Tgca observed for a reflection nebula represent a view of the dust from one direction. There is 
a very wide range of optical properties (a, g) that can flt a given observation. With one or 
more hierarchical models we can flt the best observed UV reflection nebula, NGC 7023, with 
albedos > 0.5 over the range of g that we tested (0.6 - 0.85) and with a variety of averaged 
optical depths through the model. 

In general, reflection nebulae have substantial optical depths (Text ~ 1 — 3); otherwise, 
the scattered fluxes are too faint. At substantial optical depths, uniform models can under- 
estimate the albedo rather severely because they overestimate the scattered light. Unless the 
star happens to be embedded within a dense clump, radiation can leak out of hierarchical 
clumps much more easily than from uniform dust. 

The Diffuse Galactic Light (DGL) has been interpreted (Henry 2002) with a recipe that 
predicts much more scattering from each star than our models, resulting in a low estimates 
of the UV albedo. Murthy et al. (1991) failed to detect UV DGL, but our fractal models 
often have many directions from which the star and scattered light are very faint. Thus, we 
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do not believe that the faintness of the DGL in particular directions necessarily signifies a 
low albedo. 

The determination of both a and g depends upon the variation of refiected intensity 
across the face of a centrally illuminated refiection nebula, while we have only considered 
the fiux of the scattered fight. We feel that if the fiux is as weakly constrained as our models 
show, the intensity will be similarly subject to variation because of the unknown placement 
of the star relative to the surrounding dust. 

Even the relative variation of albedo with wavelength is difficult to determine from 
reflection nebulae if there are significant changes in opacity among the wavelengths. The 
problem is that the geometry is not really the same when the opacity changes. As opacity 
per H atom increases, optically thin clumps become thick and scatter light with a different 
geometrical arrangement, so the geometry of the nebula depends on wavelength. 

Our overall assessment is that the optical properties of grains are probably as well 
constrained by theory as by observations. This statement is made in spite of well-known 
uncertainties in the theory of interstellar grains. These uncertainties are major: whether 
typical large grains are chemically homogeneous (e.g., silicate or carbonaceous) or composite, 
whether grains contain voids, or if grains have very loose ("fractal") structure. However, 
theory does not permit the large range of possible albedos provided by refiection nebulae. 
At first glance, various theories seem quite different. Weingartner & Draine (2001) have 
PAH molecules, small grains, and large grains of sificate or graphite. Li & Greenberg (1997) 
have PAHs, small grains, and large grains with silicate cores and organic refractory mantles. 
Mathis (1996) has composite grains with silicates, carbon, and vacuum intermixed within 
the same grain, plus small graphite that mimics PAHs. Draine & Lee (1984) have silicate and 
graphite with a truncated power-law distribution of sizes. Each size distribution is different. 
They all predict roughly the same a(A) because the response of the grains to radiation 
is set by fitting the well-observed interstellar extinction law Text (A), and cosmic abundances 
(though somewhat controversial themselves) constrain the materials. For instance, the values 
predicted for a(1430 A) are (0.40, 0.30, 0.37, 0.40) for Weingartner & Draine (2001), Li & 
Greenberg (1997), Mathis (1996), and Draine & Lee (1984), respectively. The A{V) are 
(0.65, 0.63, 0.53, 0.57). 

Fortunately, optical depths and column densities derived from comparisons of absorption 
and scattering along the same very narrow beam (i.e., towards a star) are not affected by 
the microstructure within the ISM. These measurements include most determinations of 

depletions of various ions in the ISM, as well as the wavelength dependence of the interstellar 
extinction law. Other implications for the hierarchical nature of dust clumping include 
a greatly increased penetration of radiation into molecular clouds and PhotoDissociation 
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Regions (PDRs). 
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Fig. 1. — The optical depth of scattering in a spherical cloud illuminated by a central star, 
plotted against the effective optical depth for the stellar radiation. Here Tgca = — ln[-Fsca/-^*o], 
Fsca is the scattered flux from the nebula, and F*o is the stellar flux if there were no reflection 
nebula. Similarly, Tcxt= — ln[-F*/F*o]. The small crosses represent the fluxes from a hierar- 
chical model with albedo, a, =0.6, phase parameter, =0.6, average radial optical depth 
for the sphere, tq, = 2, and clumping with a "fractal dimension", = 2.6. Each small cross 
represents a direction of viewing the sphere. The solid square represents ((rcxt), (Tsca)), the 
fluxes (not optical depths) of the model, averaged over viewing angles (see eqn. 1). We see 
a very wide range of points for the hierarchical model as viewed for various directions. The 
open triangle is the point for a uniform model. The circled asterisk is the observed point for 
NGC 7023, a well-observed reflection nebula, at 1300 A (Witt et al. 1993) with 40% of the 
observed extinction towards the central star to be in NGC 7023, and the circled dot is if 50% 
of the extinction is in NGC 7023. The solid lines give the relation for a uniform distribution, 
with various dust albedos marked. The dotted lines show the relationship assumed by Henry 
(2002) in analyzing the Diffuse Galactic Light. The albedos are those indicated at the nearby 
solid hues. 
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Fig. 2. — The same as in Figure 1, except that 33% of the dust density is assumed to a 
uniform density underlying the hierarchical clumps. The Tgca points are lower than in Figure 
1, showing more scattered light. 
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Fig. 3. — A plot with the same axes as Figures 1 and 2 showing the general areas of the 
(Text,Tsca) plane that are covered by points from various hierarchical models, (a, g, D, tq) 
values are: Solid line, (0.8, 0.6, 2.3, 0.5). Short-dashed line, (0.8, 0.6,2.3,0.5), with 33% of 
dust mass in constant density. Dot-dashed line, (0.6,0.6,2.6,2). Dotted line, (0.6,0.6,2.6,2), 
with constant density component. Long dashed line, (0.8,0.95, 2.3, 4). Marked points: 
observations for NGC 7023 (see Figures 1 and 2). We see that a wide variety of hierarchical 
models can cover the observations. 
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Fig. 4. — The optical depths for scattering and stellar radiation for various hierarchical 
models, averaged over viewing angles. All have (a, g, Tq) = (0.6, 0.6, 2). Each square 
represents a different arrangement of dust clumps. Panel (a) shows models with no uniform 
component; b, models with 33% of the dust in a uniform component. Filled squares: D = 
2.6; open squares, D = 2.3. Lines: the relationship for spatially uniform dust, with albedo 
as marked in panel (a), g ~ 0.6. Open circles in panel (a): dust clumped in either of two 
densities, with parameters similar to those in Witt & Gordon (2000), also with a — 0.6, g 
— 0.6, and different random spatial arrangements. The extremes of many different spatial 
arrangements are shown, along with a typical case. The dashed lines show uniform models 
with a — 0.5 and g — 0.4 instead of 0.6. We see that the other parameters are more important 
than g. 
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Fig. 5. — The squares are the same as for Figure 4, 33% of the density in underlying constant 
density, except that the averaged radial optical depth = 4 instead of 2. Open squares: D 
= 2.3. Filled squares: D = 2.6. Open triangles: purely hierarchical models (no uniform 
density), D =2.3. Filled triangles: purely hierarchical models, D = 2.6. The solid lines are 
for spatially uniform dust with a as marked, g = 0.6. 
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Fig. 6. — The albedo derived from uniform models, using the scattered flux and extinction 
of a typical purely hierarchical model as seen from various viewing angles, plotted against 
the stellar extinction. Parameters are given in the text. The "true" albedo used in the 
generating model, 0.6, is marked by the dashed line. 
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Fig. 7. — A test of whether relative values of the albedo can be derived from observations 
of a nebula at two different wavelengths or optical depths, with a given arrangement of 
dust relative to the central star. Here we compare averaged optical depths, tq, of 2 and 
1, for a particular hierarchical model. Since only tq is varied, the geometry of the dust is 
the same. The scattering and stellar fluxes at various viewing angles were interpreted with 
uniform density models. The ratio of the derived albedos is plotted against the smaller stellar 
extinction. The dashed line gives the "true" value of unity, since all points were derived with 
a = 0.6. 



